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Checking and Interpolation of Functions Tabulated 
at Certain Irregular Logarithmic Intervals 

Herbert E. Salzer 

Many functions f(x) behave as polynomials in log x. When tabulated for arguments 
in geometric progression, f(x) can be checked by ordinary differencing, and interpolation can 
be performed to a fine extent with existing tables of Lagrangian coefficients. But in practice, 
fix) is often known or calculated at some or all of the points 1, 2, 5, 10, 20, 50, 100, 200, 500* 
and 1000 (same theory for the points 0.001, 0.002, 0.005, 0.01, etc., or 0.01, 0.02,' 0.05, O.l' 
etc., or any constant multiple of 1, 2. 5, . . .). 

The present tables have a twofold use: I. Checking the correctness of fix) when tabu- 
lated at some of the more frequently occurring combinations of points 1, 2, 5, etc. This also 
includes their use to estimate the least number of tabular entries for interpolation of given 
accuracy. II. Facilitation of Lagrangian interpolation by a generalization of a scheme 
given by W. J. Taylor. 



I. Introduction 

Functions which behave like polynomials in log x 
are encountered in numerous fields, such as statistics, 
actuarial studies, economics, biometrics, electronics, 
nuclear physics, biophysics, physical chemistry, etc. 
When a function obtained experimentally is suspected 
to have this form, it is usual to examine this by plot- 
ting the values on semilogarithmic graph paper and 
observing whether they lie on a smooth curve. When 
it has been decided that the function is approximately 
a polynomial in log x (or when it is known to be ex- 
actly a polynomial in log x), the question arises as to 
how this information can be used to facilitate check- 
ing the table and interpolating in it numerically. 
This is the question to be discussed here. 

II. Arguments in Geometric Progression 

If the tabular arguments are in geometric progres- 
sion, such as 1, 2, 4, 8, . . .; or 1, 10, 100, 1000, . . „ 
the problems of checking and interpolation are quite 
simple, because the function is effectively tabulated at 
a constant interval in log x. Thus an examination of 
the ordinary differences of the tabular values will 
reveal any errors, as well as the lowest degree which 
an approximating polynomial in log x must have in 
order to yield a certain desired accuracy. Either the 
differences themselves can be used for interpolation, 
or if we wish, to avoid using differences, any one of 
several well-known tables of Lagrangian interpola- 
tion coefficients may be used for interpolation to a 
very fine subdivision of the tabular interval. The 
most extensive of these tables of interpolation co- 
efficients are contained in [I]. 1 

III. Divided Differences for Checking 

In practice, however, it is often found that the 
function has been determined at some or all of the 
points 1, 2, 5, 10, 20, 50, 100, . . ., 1000 (or a con- 
stant multiple of those values, e, g., 0.1, 0.2, 0.5, 
1.0, . . .). Then it is no longer possible to use the 

1 Figures in brackets indicate the literature references at the end of this paper. 



ordinary differences for checking purposes. Instead 
we use a certain generalization, known as "divided 
differences/' a subject that is treated fully in most 
textbooks on finite differences (see [2] or [3]). The 
(n— l)th divided difference with respect to log x, of 
a function f(x) tabulated at the n points x u x 2 , . . ., 

n 

x n , can be written in the form y^A/fe), where the 

A t are certain numbers depending only on the x t , 
not on the function. It is convenient^ to tabulate 
these quantities, once and for all, for the usual sets 
of values of x t and for the various values of n, 
so that the divided difference can be obtained by 
a single accumulation on a calculating machine. 
This has been done here for 72=3(1)10 and x t ranging 
from 1,2, 5, . . ., to 1000, in table 1, A; for n=3 (1)7 
and xt ranging from 1, 5, 10, . . . to 1000 in table 
1, B; and for ^=3(1)7 and x t ranging from 1, 2, 10, 
. . . to 1000 in table 1, C. All entries are given to 
eight significant figures, correct to about a unit in 
the last place. 

The explicit expression for A t is 



:1/V QogioZ* — logioZy), 

i 



(1) 



where w' indicates the product over alljVi. It can 

3 

be shown that if/(z) is a polynomial of degree n— 2 in 
log x } then the {n— l)th divided difference is identi- 
cally zero. If the (n— l)th divided difference is 
sufficiently small, then/(;r) considered as a polyno- 
mial of degree n~ 2 in log z, is in all likelihood, free 
from error. (However the user is cautioned that a 

n 

sufficiently small value of ^ AtfCxt) is only a strong 

i=l 

indication, rather than conclusive evidence, of the 
correctness of the n entries, considered as values of a 
polynomial of degree n—2.) Now the tabular entries 
may be entirely correct and the user might still want 
to employ the most economical Lagrangian formula. 
Again, this information is conveyed in the sufficiently 
small (n— l)th divided difference, which establishes 
the adequacy of an (n— 1) point Lagrangian inter- 
polation formula. In brief, if the function is checked 
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Table 1. Coefficients A t for checking and interpolation, based upon the points Xi 

A 



x% 


Ai 


Throe-point 


Four-point 


Five-point 


Six-point 


X\ = l 
2-2 = 2 
2-3 = 5 
J4 = 10 
2:5 = 20 
2*6 = 50 


Ax = 
A 2 = 
A%— 
Ai = 

Ah = 

A 6 = 


4. 75260 47 

-8. 34781 13 

3. 59520 66 


-4. 75260 47 

11.94301 8 

-11.943018 

4. 75260 47 


3. 65295 55 
-11.94301 8 

19. 83692 3 
-15.78781 1 

4. 24095 00 


-2.15010 00 

8. 54329 79 

-19.83692 3 

22. 58725 1 

-10.65726 

1. 51373 42 










A— continued 


Xi 


Ai 


Seven-point 


Eight-point 


Nine-point 


Ten- point 


2-, = l 

2-2 = 2 
2:3 = 5 
2-4 = 10 
2-5 = 20 

J6 = 50 
J7-100 

2-8 = 200 

2:9=500 

r 10 — 1000 


.li = 
A 2 = 
Az = 
Ai = 
A,= 

,16 = 

-1: = 

.4 8 = 
yU = 
A io = 


1. 07505 00 
-5. 02851 60 
15. 24709 1 
-22. 58725 1 
15. 24709 1 

-5.02851 60 
1. 07505 00 


-0. 46720 383 
2. 51425 80 
-9. 51717 89 
17. 36105 3 
-15.24709 1 

8. 35218 43 

-3. 57123 89 

0. 57521 751 


0. 17310 449 
-1.04850 75 

4. 75858 94 
-10.21857 5 

10. 90682 8 

-H. 35218 43 

5. KI928 77 
-1.44548 80 

0.11694 543 


-0. 05770 1497 
0. 38848 430 

-2.06802 58 
5. 10928 77 

-6. 41967 08 

6.41967 08 

-5. 10928 77 
2.06802 58 

-0. 38848 430 
0.05770 L497 



















Xi 


.1, 


Three-point 


Four-point 


Five-point 


Six point 


Seven-point 


X\ = 1 

.T2 = 5 
23 = 10 
X4 = 50 

2' 5 =100 
J6 = 500 
2-7 =H)(!0 


/ll = 
.12 = 

.4 3 = 

.4 4 = 

.45 = 

.4 6 = 
-1: = 


1.43067 66 

-4.752(11 17 
3.32192 81 


-C. 84208 465 
4. 75260 47 

-1.75260 47 
0.8420S 465 


0. 42104 232 
-3. 65295 55 

4. 75260 47 
-2. 79734 47 

1.27665 32 


-0. 15600 111 
L. 82647 77 

-2.79734 47 
2. 79734 47 

-1.82647 77 
0. 15600 111 


0. 05200 0371 
-0. 79376 529 

1.39867 23 
-2. 15010 00 

1.82647 77 
-0. 51822 448 

0. 18493 938 























Xi 


.1, 


Three-point 


Four-point 


Five-point 


Six-point 


Seven-point 


2:1 = 1 

2:2=2 

2:3=10 

2:4=20 

2:5=100 

2:6=200 

2-7=1000 


.4i = 
.4 2 = 

.4 3 = 

-44 = 

-4 5 = 
.4 6 = 
A 7 = 


3. 32192 81 

-4. 75260 47 

1.43067 66 


-2.55330 63 
4. 75260 47 

-4. 75260 47 
2. 55330 63 


1.27665 32 
-2.79731 47 

4. 75200 47 
-3. 65295 55 

0.42104 232 


-0. 55481 813 
1.39867 23 

-3.65295 55 
3.65295 55 

-1.39867 23 
0. 55481 813 


0. 18493 938 

-0. 51822 448 
1.82647 77 

-2. 15010 00 
1. 39867 23 

-0. 79376 529 
0. 05200 0371 
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by the A/s for a certain n, we use a different set of 
A/s corresponding to n— 1, for interpolation. 

IV. Interpolation 

Once we have ascertained, by use of table 1, that 
a function is adequately represented by a polynomial 
of degree n— 1 in log x (remembering that n— 1 is 
here the /? — 2 of III.), we can interpolate by use of a 
rearrangement of Lagrange's formula, which was 
suggested by a paper of W. J. Taylor [4]. The original 
method and notation of Taylor are not described in 
this article because he developed them only for the 
special case of equally spaced arguments, and here 
they would be superfluous. The interpolation 
formula is 



where 



/W-(g«</W)/g^ (2) 

ai = ^/(log x — log x^ (3) 



and A t is given by (1). 

V. Logarithms to Other Bases 

The coefficients A t have been computed from 
common (base 10) logarithms, according to (1). It 
would have been possible, instead, to have the A/s 
calculated for any other logarithmic base. Further- 
more, it is permissible when computing the a/s from 
(3), to use logarithms to a base different from that 
underlying the A/s. To illustrate this point, natural 
logarithms (to the base 6 = 2.71828 183) have been 
used in the two examples in section VIII. Extensive 
tables of log e £ are given in [5]. 

VI. Conversion to Arguments in Geometric 
Progression 

If a large number of interpolations are required 
for a function given for any one set of values x u it 
may be convenient to prepare, by the method in 
section IV, an auxiliary table giving the function at 
a new set of values y u where the y f are now in geo- 
metric progression. Then we can use ordinary in- 
terpolation formulae and coefficients, as suggested 
in section II. This point was called to the writer's 
attention by Churchill Eisenhart and Julius Lieblein, 
of the Statistical Engineering Laboratory. 

VII. Use of Ai for Other Arguments 

If the given set of values x t does not begin with 1, 
the table of coefficients A t may sometimes still be 
used by a suitable change of the independent varia- 
ble x. Notice that a polynomial in log (ax) or log 
Q>jx) is still a polynomial in log x. Thus, for instance, 
if a function is given at the points #=10, 20, 50 and 
100, we may consider it as a function of x' = 100/x, 
given at the points x f = \,2, 5 and 10, so that table 1,A, 
may be used. Again, if a function is given at the 
points #=20, 100, 200 and 1000, it may, instead, be 



considered as a function of x r =1000/# given at the 
points x' = l, 5, 10 and 50, so that table 1,B, may be 
used. Schedules A, B, and C list such transforma- 
tions of the independent variable, giving both the 
sets of x t for which the change of variable to 1/x is 
applicable, and the new arguments #/, which are 
proportional to l/z*. 

Transformation schedules 
Schedule A. 



Xi Xi 


Xi 


Xi 


Xi 


Xi' 


X{ 


Xi 


Xi 


Xi' 


2 


5 


2 


50 


2 


500 


5 


20 


5 


200 


5 


2 


5 


20 


5 


200 


10 


10 


10 


100 


10 


1 


10 


10 


10 


100 


20 


5 


20 


50 






20 


5 


20 


50 


50 


2 


50 


20 






50 


2 


50 


20 


100 


1 


100 


10 






100 


1 


100 

200 

500 

1000 


10 
5 
2 
1 






200 

500 

1000 


5 
2 
1 



Schedule B. 



Xi 


Xi' 


Xi 


Xi 


Xi 


Xi' 


Xi 


x/ 


2 


50 


2 


500 


5 


10 


5 


100 


10 


10 


10 


100 


10 


5 


10 


50 


20 


5 


20 


50 


50 


1 


50 


10 


100 


1 


100 


10 






100 


5 






200 


5 






500 


1 






1000 


1 











Schedule C. 



Xi 


Xi 


Xi 


Xi 


Xi 


Xi 


Xi 


Xi' 


5 


20 


5 


200 


2 


10 


2 


100 


10 


10 


10 


100 


10 


2 


10 


20 


50 


2 


50 


20 


20 


1 


20 


10 


100 


1 


100 


10 






100 


2 






500 


2 






200 


1 






1000 


1 











VIII. Illustrations of Use of Tables 

The following two illustrative examples show how 
to use the tables of A t . 

Example 1. Given 



Xi 


/(*<) 


Xi 


f(Xi) 


1 


0. 52 


50 


6333. 56 


2 


1. 11 


100 


15752. 58 


5 


50. 30 


200 


34648. 66 


10 


343. 74 


500 


85417. 78 


20 


1447. 15 


1000 


155833. 72 
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To calculate j (18). In the absence of knowledge of 

an explicit formula for f(x), we find that its seventh 
divided difference, employing the coefficients A t in 
table 1, A, for the eight points from 1 to 200, is less 
than 0.05 in absolute value. This indicates that the 
function behaves, up to the last place, like a sixth 



degree polynomial in log r (any base). Hence only 
seven points are needed for the interpolation. The 
work may be arranged in the following way: (Note 
that an extra place is carried in some columns to 
reduce the error that would arise in rounding before 
the final answer is reached.) 



Xi 


log c 18 — log e Zi 


a;= AJ (log e lS — logeZ;) 


Oiffai) 


1 

2 

5 

10 

20 

50 

100 


2. 89037 176 

2. 19722 458 

1. 28093 385 

0. 58778 667 

--0. 10536 051 

-1. 02165 125 

-1. 71479 843 

* 


0. 371942 

-2. 288576 

11. 903106 

-38. 427634 

-144. 713527 

4. 921950 

-0. 626925 


0. 193 

-2. 540 

598. 726 

-13209. 115 

-209422. 181 

31173. 466 

-9875. 686 


7 

]T>.-=-168. 859664 
i = i 


7 

Sa./W = - 200737. 137 
i =i 








S a »-/fo)/S°»-= 1188 - 78 



The answer found here, 1188.78, is correct to its 
last significant figure. The function /'(.r) was chosen 

to be (log e z) 6 +3(log € £) 5 +f 



Exam/pie 2. Given 



6 



Xi 


Xi' 


f(Xi) or./V) 


20 

100 

200 

1000 


50 

10 

5 

1 


15. 45981 
40. 07131 
54. 18217 
95. 099 19 



To calculate J (160). To use the present tables of 
coefficients A* it is necessary to change the variable 
to z' = 1000/.r. (See Transformation schedule B.) 
If /(z) is called j(x') 7 then /(160)==/(6.25). As in 
the previous example, in the absence of knowledge of 
an explicit formula for /(./•), we find that the absolute 
value of its third divided difference, employing the 
coefficients A, in table 1, B, for four points, is less 
than 1% units in the fifth decimal place. This indi- 
cates that the function behaves, up to the last place, 
like a second degree polynomial in log x. Hence 
only three points are needed for the interpolation, 
which is carried out as before: 



Xi' 


loga 6.25- log e Xi' 


en :=Ai/ (loge 6.25 -log, x^) 


ajfa') 


1 

5 

10 


1. 83258 146 
0. 22314 355 
0. 17000 363 


0. 78068 92 

-21. 29841 84 

-7. 06787 75 


74. 24314 

-1153. 99453 

-283. 21911 


][>i=-27. 58560 67 
i =i 


3 

XX / (»<') = - 1362. 97050 
i = i 

2>*7V)/s o <= 49 - 4 ° 876 



This answer, 49.40876, is also correct to the last 
figure given, the function f(x) having been chosen 
to be 2.20(log e x) 2 — 1.43 \og e x. 



The author expresses his appreciation to Churchill 
Eisenhart and Julius Lieblein for their thorough 
study and constructive criticism of this paper in its 
original form, and to John Todd and Franz L. Alt 
for their editorial review and for suggesting some 
final improvements in the presentation. 
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